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Abstract. The evolution of the radio emission of shell-type Supernova remnants (SNRs) is modeled within the 
framework of the simple and commonly used assumptions that the mechanism of diffusive shock acceleration 
(DSA) is responsible for generating radio emitting electrons and that the magnetic field is the typical interstellar 
field compressed at the shock. It is considered that electrons are injected into the mechanism in test-particle 
regime directly from the high energy tail of the downstream Maxwellian distribution function. The model can be 
apphed to most of the observed SNRs because the majority of detected SNRs are shell-types and have a more or 
less spherical shape and are sources of nonthermal radio emission. It is shown that the model successfully explains 
the many averaged observational properties of evolved shell-type SNRs. In particular, the radio surface brightness 
^ ' (E) evolves with diameter as ~ £)-(o.3-o.5) ^ while the bounding shock is strong (Mach number is 7W > 10), 

^\ I followed by steep decrease (steeper than ~ _D~'''^) for M < 10. Such evolution of the surface brightness with 

. diameter and its strong dependence on the environmental parameters strongly reduce the usefulness of E — D 

' relations as a tool for determining the distances to SNRs. The model predicts no radio emission from SNRs in the 

QQ I late radiative stage of evolution and the existence of radio-quiet but relatively active SNRs is possible. Our model 

easily explains very large-diameter radio sources such as the Galactic Loops and the candidates for Hypernova 
radio remnants. The model predicts that most of the observed SNRs with Eighz m ^ sr ^ Hz ^ are 

located in a tenuous phase of the ISM. The model also predicts the existence of a population of 150 — 250 pc SNRs 
with EiGHz ^ lO"'^^ Wm"'^ sr~^ Hz~^ if the kinetic energy of the explosion is ~ lO^'^ erg. From the comparison of 
Oh' the model results with the statistics of evolved shell-type SNRs, we were able to estimate the fraction of electrons 

accelerated from the thermal pool in the range (3 -^ 11) x 10~*. If acceleration takes place directly from the high 
^ • energy tail of the downstream Maxwellian distribution function, then the corresponding injection momentum is 

^ I estimated as pinj ~ (2.7 — 3) • pth 
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^ 1. Introduction dio emission from the SNR shells has still not been fully 
- - ' solved. Since the pioneering works on particle acceleration 
The radio synchrotron radiation is a prominent output colHsionless shocks by Krymskii (1977), Bell (1978a), 
of SNRs and is normally used to identify them among ^^^^ Blandford & Ostriker (1978), the DSA mechanism 
other radio sources. Although SNRs are studied at almost ^^^^ ^^^^^ regarded as the most suitable mechanism re- 
every wavelength ranging from gamma-rays through low- sponsible for the origin of radio-emitting electrons in shell- 
frequency radio waves, the radio observations still remain Although we directly observe only electrons, 
the most effective means of obtaining important informa- mechanism of DSA has theoretically more successful 
tion on these objects; m fact, most of the Galactic SNRs application to the protons and heavy ions. The effective- 
have been discovered by their radio emission. Naturally, ^^^^ acceleration for ions may be very high in certain 
much richer observational information on the SNRs is parameters regimes (see the review by Malkov & Drury 
available m the radio range; m the last version ( January 200I), and this mechanism is believed to be the source of 
2004) of the Green (2004a) catalogue 231 remnants are g^j^^^j^ ^^^^-^ ^^^^ ^-^^ ^^^^^-^^ ^p ^qis^v. In the 
hsted, a number that is continuously growing (e.g., Brogan ^^^^ electrons, not only is the electron acceleration the- 
et al. 2006). For most recent review of the radio SNRs, see o^etically problematical but the mechanism for electron 
Reich (2002). heating in coUisionless shocks is also poorly understood. 

The question of the origin of the high-energy electrons ^t the same time, synchrotronically radio-emitting rela- 

and magnetic fields responsible for the synchrotron ra- tivistic electrons in shell-type SNRs with spectra coinci- 
dent with the prediction of the mechanism can be con- 
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sidered as the only direct evidence of a DSA action. The 
recently observed X-ray synchrotron emission from several 
SNRs has strengthened this evidence (see Ballet 2005 for 
a recent review). 

There are a number of models describing the evolution 
of the radio emission of yoimg SNRs (e.g., Chevalier 1982; 
Gull 1973; Dickel et al. 1989) and of the SNRs in cloudy 
environments (e.g., Chevalier 1999; Bykov et al. 2000). 
The most popular model for an interpretation of the ra- 
dio emission from evolved and old SNR is the mechanism 
of van der Laan (1962) and its modifications. According to 
this model the magnetic fields and the interstellar cosmic 
ray electrons compressed behind a cooling shock front are 
responsible for the continuum synchrotron radio emission 
of SNRs, so this mechanism is applicable to the remnants 
in the radiative phase of evolutions. In spite of the fact 
that the general theory of SNR predicts the inevitable 
onset of a radiative phase in the course of SNR evolution 
(e.g., Chevalier 1974), observationally very few candidates 
for radiative SNRs are detected. The set of SNRs, which 
have long been considered as highly evolved remnants in 
the radiative phase, include IC 443, W44, CTB 80, S147, 
MSH 11-61A, OA 184, etc. However, the intensive op- 
tical emission observed from these SNRs in most cases 
can be interpreted by the interaction of a non-radiative 
adiabatic shock wave with the dense interstellar clouds. 
The co-existence of intensive optical emission and high- 
temperature plasma, which gives the thermal X-rays, is 
considered as an argument in favor of such interpretation. 
Indeed, in several SNRs there is direct evidence of slow 
radiative shocks, which are driven into ambient clouds by 
a fast nonradiative blast wave (McKee & Cowie 1975). For 
instance, in Seta ct al. (1998) clear evidence of the inter- 
action of SNRs W44 and IC 443 with the molecular clouds 
is presented. In the Cygnus Loop SNR, the bright optical 
filaments are attributed to shocked clouds (e.g., Hester & 
Cox 1986; Patnaude et al. 2001). Shocked clouds have also 
been observed in Vela SNR (e.g., Bocchino ct al. 1999), 
so there are very few candidates for SNRs that could be 
considered as real radiative remnants in accordance with 
the standard theory. The S147 SNR (Reich 2002) can be 
considered as a representative of this tiny group of SNRs. 

The lack of radiative SNRs in our galaxy can be ex- 
plained in several ways. First, there are inevitable selection 
effects, mainly due to the irregular high-intensity back- 
ground emission, which works against the detection of ex- 
tended SNRs evolving in the rarefied environs. Second, as 
shown (Blondin et al. 1998; Bertschinger 1986; Vishniac 
1983; Falle 1981), the shells of SNRs at the beginning of 
the radiative stage become highly unstable leading to frag- 
mentation of the dense shell, thereby making the action of 
the Laan's mechanism become inefficient. And, finally, if 
the density of the medium in which the SNR is located is 
low enough, it is possible for the SNR to finish its life by 
merging with the ISM before cooling becomes important. 
In addition, it is important to take into account the ef- 
fect of accelerated particles, trapped in the matter behind 
the shock front, which may impede the processes of catas- 



trophic cooling. The problem of the visibility of old ra- 
diative SNRs in the Hi 21-cm emission line in our Galaxy 

was considered in the recent work by Koo & Kang (2004) , 
who also noted a considerable deficit of radiative SNRs in 
our Galaxy. 

As a result, the van der Laan mechanism encounters 
a number of difficulties in interpreting the radio emission 
from the observed large-diameter SNRs. However, under 
some specific conditions, the modification of this mech- 
anism proposed by Blandford & Cowie (1982) (see also 
Bychkov 1978) can be effective in generating radio emis- 
sion from old and large SNRs. Although these models also 
include the mechanism of DSA, the specific cloudy ambi- 
ent interstellar medium is required to be effective in gener- 
ating radio emission; besides, the physics of cloud - shock 
interaction is far from being understood in detail. 

Though the ISM generally has an inhomogeneous na- 
ture, many extended and old SNRs, as a rule, demonstrate 
nearly circular shells in radio and x-ray, which implies 
that they have been expanding into a relatively homo- 
geneous region of the ISM. In practice, the circular shape 
is one of the observational properties of SNR by which 
they are identified as SNR. In such homogeneous regions, 
it is possible, however, that such small-scale clouds ex- 
ist (Nagashima et al. 2006), as required by the model 
of Blandford & Cowie (1982). Two cases are possible. 
If small-scale clouds have high spatial concentration, the 
evolution of SNR occurs according to the scheme con- 
sidered by Cowie et al. (1981), which will result in the 
morphology being completely different from the standard 
Sedov solution (see also Chieze & Lazareff 1981). Indeed, 
a small fraction of SNRs in X-rays show the morphol- 
ogy consistent with such a possibility. According to Rho 
& Petre (1998), up to 25% of all X-ray-detected Galactic 
SNRs show this "mixed-morphology" : i.e., center-filled X- 
ray and shell-like radio morphology. Moreover, in such a 
medium the SNR cannot extend up to the large diameters, 
which allows us to exclude this case from consideration. 
In the other case of a low concentration of the interstellar 
cloudlets, their contribution to the integral radio emis- 
sion will be minor, and the remnant will have a more or 
less regular shell morphology. In this case the mechanism 
responsible for generating the radio-emitting relativistic 
electrons works at the outer shock front that regularly 
bounds the SNR. In reality, we do not exclude the possi- 
bility that the van der Laan mechanism (1962) will work 
locally in some parts of the large diameter remnant; in 
other words, different stages can occur simultaneously in 
different locations within a single remnant. Very proba- 
bly this is the case in above-mentioned, optically bright 
galactic SNRs. 

Recently Berezhko & Volk (2004) developed a time- 
dependent model of radioemission from shell-type SNRs, 
taking into account the non-linear shock acceleration ef- 
fect and the generation of the magnetic field in the shock 
vicinity. This model is very good for young SNRs, but 
in case of evolved and old remnants, the application of 
both of these effects can be questioned. Therefore, it is im- 
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portant to construct a model that is based on the very 
general results of the DSA theory, and that resorts to 
minimal non-standard assumptions. In this paper wc sug- 
gest the model describing the evolution of radio emission 
of SNRs, which is based on the very common assump- 
tions that the radio-emitting electrons are accelerated dif- 
fusively at the shock front from the thermal pool in the 
test-particle regime. The aim of the work is not to model 
a specific SNR in detail but to illustrate the general be- 
havior of the radio properties of SNRs (surface brightness, 
spectral index, etc.) with their evolution, to test the capa- 
bility of DSA in reproducing the general statistics of the 
shell-type SNRs, and to constrain the injection problem 
as far as possible. 

The set of main assumptions are the following: 

- We consider spherical adiabatic SNR evolving in the 
homogeneous ISM. To follow the evolution of the remnant 
up to the maximum sizes, we used an analytical approx- 
imation of Cox & Andersen(1982), which describes the 
evolution of an adiabatic, spherical blast wave in a ho- 
mogeneous ambient medium of finite pressure. At early 
times while the shock wave is strong, this approximation 
coincides well with the zero-pressure self-similar Scdov so- 
lution. In applications of the model to real SNRs, the pres- 
sure of the interstellar magnetic field is included. 

- Electrons are accelerated diffusively at the main 
shock front, and their injection takes place from the down- 
stream thermal distribution function. The threshold en- 
ergy of injection is proportional to the temperature of the 
downstream plasma, and as a conseqiicnce of this assump- 
tion, the injection of new particles in the process of accel- 
eration stops when radiative cooling of the shocked gas 
first starts. The test particle approximation is used, all 
the nonlinear effects are neglected, and electron-ion tem- 
perature equilibrium is presumed. It is also assumed that 
accelerated electrons only lose their energy adiabatically 
and there is no diffusion. 

- The magnetic field is the typical interstellar field 
compressed at the shock wave up to four times; i.e. there 
is no additional amplification of the interstellar magnetic 
field in the shell of the SNR. It is also assumed that it 
is fully chaotic and is frozen-in to the plasma inside the 
remnant. 

The paper is organized as follows. In Sect. 2 the 
radially-dependent distribution function of accelerated 
electrons is derived, which was then used to obtain the ba- 
sic formulas describing the radio properties of SNR. The 
list of main predictions of the model is presented at the 
end of the section. To apply our model to real SNRs the 
initial conditions and parameters are specified in Sect. 3. 
Further, the empirical E — £> relation is constructed to 
compare the predictions of the model with observational 
data . With the help of this relation, the values of the in- 
jection parameter is estimated. The final section contains 
our main conclusions. 



2. Model 

2.1. Distribution function of accelerated electrons 

In the test particle approximation, the action of DSA in 
planar, steady shock front leads to an isotropic momentum 
distribution function for the accelerated particles of the 
form 



/a = « [H{j> - Pinj) - H{p - Pmax)] 



(1) 



where A is the normalization constant, H the Heaviside 

step function, pinj the injection momentum, and Pmax 
the maximum momentum of electrons. A very important 
property of this spectrum is that the spectral index q, de- 
fined as 



(2) 



depends only on the shock compression ratio Xs = p/ po, 
where po and p are mass densities for the up- and down- 
stream states, respectively. As a parameter describing the 
evolution of SNR in the following we mainly use the Mach 
number defined as the ratio of shock velocity Vs to the 
sound speed (correctly, to the maximum signal velocity 
ahead of the shock) in the ISM Cg, M.^ = Vg/Cg. From 
the standard theory of the gas with an adiabatic index 
7 = 5/3, the compression ratio at the shock is expressed 
as (Landau & Lifshic 1986) 



4Mi 



(3) 



and for the dependence of the spectral index q on the Mach 
number, we have 



(4) 



It is important to note that it is much easier to apply 
the mechanism of DSA to the evolved SNRs than to the 
young ones. This is first because it is expected that various 
nonlinear and geometrical effects do not play an essential 
role in the evolved SNRs, so they can be neglected and 
the test particle approximation is acceptable. Second, the 
structure of the shock wave of the evolved SNRs in homo- 
geneous ISM can be described in a simple way, unlike the 
young ones. In the young SNRs the shape of the circum- 
stellar matter distribution often plays a crucial role, and 
they often have a double- or even multiple- shell struc- 
ture. Third, it is expected that in the evolved SNRs the 
temperature equipartition between electronic and proton 
components takes place. 

One very important problem remaining still unsolved 
in the thcjory of DSA is the injection of particles into the 
action of the mechanism. Although in the case of pro- 
tons, there is some progress in this direction (e.g., Drury 
& Malkov 2001), the theory of injection practically is not 
advanced for electrons. One basic theoretical difficulty in 
applying DSA to electrons is that it is still unclear how 
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thermal electrons gain momentum above the thermal ions, 
so they can be injected into a population that can respond 
to DSA. In our case, the determination of the coefficient A 
in the spectrum (1) is directly connected to this injection 
problem. As initially done by Bell (1978b) in applications 
of DSA to real objects, the injection problem is usually 
avoided by assuming that some fraction of the thermal 
particles are "injected" as non-thermal particles at some 
"injection momentum". pi„j, which separates the thermal 
particle population from the non-thermal. Here we intro- 
duce these two parameters as 



V = Ninj /rie 

V' = PinjMh 



(5) 
(6) 



where TVjnj is the concentration of electrons contained in 
the spectrum (1), pinj the injection momentum, and pth 
the thermal momentum of heated by the shock plasma 
with number density Uc- 

In present paper we assume that suprathcrnial elec- 
trons from the tail of the downstream Maxwellian distri- 
bution function are subjected to acceleration, which allows 
us to determine the coefficient in the spectrum of the ac- 
celerated electrons through the equality 



where 



/m(p) 



exp - — 
Pth 



is the Maxwellian distribution function of electrons heated 
and compressed by the shock, and pth is the thermal mo- 
mentum of electrons pth = (2mers)^/^. We also assume 
that there is an equilibrium between electron and proton 
temperatures, i.e. Te = Tg = (3/16) furipVg {j2 is the mean 
mass per particle in units of the proton mass), which is 
expected to be the case in evolved SNRs. Prom the above 
equality, we have for A an expression 



A = 



47r(l - £9- 



,9-3 



where £=pi^Jp 

77= i = — 



(7) 



(8) 



with 77 « 4.60 X IQ-^, 5.85 x lO'^, 4.40 x lO"" for 
ip = Pinj/Pth = 2, 2.5, 3, respectively. In deriving expres- 
sion (7) we have used = Xs Ucq from the shock jump 
conditions and q/ (g — 3) = Xg, which follows from (2). 

Taking (7) into account, we can rewrite the spectrum 
of accelerated electrons as 



fa 



47rpf„j(l-e9-3)%Pi„J 

[H {p - Pinj) - H{p- Pma: 



(9) 



It is interesting to note that this distribution function has 
the same form as the one derived by using the 5 - function 
as the source of injected particles. 

This elementary "solution" of the problem of injection 
has some theoretical basis. Indeed, theoretical works de- 
voted to the problem of electron injection and following 
acceleration in various approximations show that the over- 
all spectra of electrons has modified Maxwellian distribu- 
tions with power-law tails at high energies (e.g., Bykov 
& Uvarov 1999). The hypothesis used above about the 
injection of particles directly from the Maxwellian distri- 
bution is not entirely physically correct. Only under very 
specific conditions can superthermal electrons be accel- 
erated by DSA, even though in this case our considera- 
tion is very simplified. The problem of particle accelera- 
tion from the Maxwellian distribution was considered in 
detail by Gurevich (1960) (see Bulanov & Dogiel 1979 
for the case of shock acceleration and Dogiel 2000), who 
showed that deviation from the equilibrium Maxwellian 
distribution begins at the momenta much lower than pi„j , 
and there is broad transition region between the equi- 
librium Maxwellian distribution and power-law spectrum 
at high energies. As shown by Asvarov et al. (1990), 
this part of the electron spectrum may manifest itself 
through the bremsstrahlung X-ray emission at photon 
energies ^ 10 keV. Indeed, such emission has been de- 
tected from several SNRs, although most authors inter- 
pret it as synchrotron emission from TeV-energy electrons; 
however, such an interpretation cannot be considered as 
proved, especially in the case of Gas A SNR (e.g., Vink 
2003). Physically it is most probable that the low energy 
suprathermal electrons are accelerated first by the electric 
fields (or stochastically) and then by DSA to the higher en- 
ergies. Although it is difficult to offer any concrete mecha- 
nism accelerating electrons from the thermal distribution 
up to the energies sufficient for DSA, but our simple as- 
sumption can be considered as the first approximation, 
and it is more or less good if ^ noo- And finally, if we 
accept this naive representation about injection then the 
two Bell parameters turn out to be related to each other 
through relation (8). 

The distribution function (9) is derived for the plane 
and steady shock, neither of which is fulfilled in the case of 
SNRs. To take the effects of nonstationarity and sphericity 
into consideration and to follow the evolution of the spec- 
trum of accelerated particles in a simple way, it is useful 
to employ the method known as "onion-shell" (Moraal & 
Axford 1983; Bogdan & Volk 1983). For our purposes, the 
method of Moraal & Axford (1983) is more appropriate, 
so we use it with some modifications. 

First of all, we assume that at any moment of interest 
the spectrum of electrons (9) extends up to the energies 
(say, 30 GeV) sufficient for emitting in the range of radio 
frequencies (0.1 -j- 10) GHz for the plausible values of the 
magnetic fields between 3 fiG and 10~* G, in other words, 
in formula (9) we consider that Pmax ^ oo (e = 0). 

Taking these assumptions into account and express- 
ing the injection momentum as a function of shock Mach 
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number pi„j = poi{-Ms/ Aio), where poi is the injection 
momentum at the initial value of the Mach number, A4o, 
when the shock radius was Rq, the spectrum (9) takes the 
form 



/a = Trrr 9 I — 



477 p^i 



Mo 



q-3 



H 



P- 



«,^ ) (10) 



Let us introduce the radially dependent distribution 
function, f{p, R, Rg), which will allow us to find the con- 
centration of accelerated particles at any radial position 
inside the SNR with a given radius Rg. To do this, one 
must first adopt a concrete model describing the struc- 
ture of the SNR. As was noted before for the model of 
the SNR, we use an analytical approximation of Cox & 
Andersen (1982) (details of this approximation are given 
in Appendix). Second, we assume that the accelerated 
electrons in the blast wave shell are frozen-in to the matter 
and participate in the bulk motion of the gas without any 
dynamical effect; i.e. there is no diffusion of electrons. This 
is not a very strong constraint for the value of the diffusion 
coefficient, if we take into account that in the gyroradius 
limit for ^ 1 GoV electrons in the B ~ 10~^ G magnetic 
fields Kmin — c^p/3eB ~ 3 x 10^^ cm^/s, and there is ob- 
servational evidence that the actual value of the diffusion 
coefficient in the region of the shock front is only a few 
times larger than Kmin (Achterberg et al. 1994). Finally, 
we consider adiabatic cooling of the accelerated electrons 
in the course of expansion of the shell. The number of par- 
ticles, dN, in the momentum range {p;p + dp) contained 
inside the SNR with R^ is 



Its 

dN = 4Trp^dp j 4'KR'^dRf{p,R,Rg), 



(11) 



where i?in is the new location of a mass element that was 
first shocked at Rq. We assume that there are no particles 
inside the radius Rin, which means that the action of the 
acceleration mechanism first starts when the shock radius 
is i?o- 

The same number of particles, dN, equals the sum of 
particles that were initially accelerated into the momen- 
tum range {pi,Pi + dp\) and were contained in the com- 
pressed volume dVi = AnR^^dRi/xsi {xgi is the compression 
ratio at the moment when shock radius was Ri) and sub- 
sequently cooled down to {p:p + d,p). To take the adiabatic 
cooling into account, we introduce a cooling function as 
in (Moraal & Axford 1983) 



C = p/pi = 



1 dVi 
Xsi dV 



1/3 



which describes loss of particle momentum from pi to p 
due to expansion of the mass element from 1/pi to 1/p. 
Using (A. 3) and (A. 4), the cooling function takes the form: 



C={x{r,Rs)/xgif\ 



(12) 



Since the isotropic distribution function remains un- 
changed by this type of cooling, now we can find dN by 
integrating over all volumes dVi 



dN = Awp^dp J inRfdRi 



U(p/C,Ri) 



(13) 



where C depends on Ri and Rg. Equating (11) and (13), 
we have 



r'^dr ■ f{p, r, Rg) 



Rq/ Re 



where we introduced notations for the normalized radii as 

r = R/Rs and n = Ri/R^. Using the expressions (A. 3) 
and (A. 4), after transformation r; = ri(r). lirj = = 

2 

x{r) ^dr in the right hand side integral and taking (10) 
into account, we obtain 



/(p,r,ii,) = ^giC* ^ 

47rp3i \po 



P 



Mo, 



H{r-nn)H [p-Cp, 



Mi 
Mo 



(14) 



In this spectrum, gi = qi{r,Mi) is the spectral index of 
particles that were accelerated earlier when the shock front 
was at radius Ri with strength Mi and that are now at 
R = r Rg inside the SNR with radius Rg. For any radial 
position r at the present shock strength Mg, the value of 
Mi can be found by solving the equation (see Appendix 
1) 



{Mf - 11/15) 
{Mf 1) 



2/3 



ri{r,Mg) 



{Ml ~ 11/15) 
(M^s - 1) 



2/3 



(15) 



The value of inner radius rin in (14) for any Rg (and Mg, 
respectively) can be found as the root of Eq. (15) but with 

substitutions of Mo for Mi and Hn for r. 

Now if we write the distribution function in the form 



dN = K^E^^'^-^^ dEi, (electron/cm^), 



(16) 



than taking into account an equality Kf.E^''''~'^^ dE\. = 
Anp^dp f{p,r, Rg) and that E = pc for the range of en- 
ergies of interest for the coefficient in (16), we obtain an 
expression 



jre = ?7neo(poic)*-'«C* 



Mi 
Mo 



9i-3 



X H{r-ri^)H{p-CpoiMi/Mo). 



2.2. The coefficient of synchrotron emission. 



(17) 



Isotropically distributed electrons with the distribution 
function (16) in the fully chaotic magnetic field B radiate 
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synchrotron radio emission with the coefScient (Ginzburg 
1981) 



and 



£^ = 1.35 X 10-22 a(a) ifeB 



1 fQ.26 X lO^S' 

3_-lxj„-l> 



(erg cm" sr"^ Hz"') (18) 

where v is the frequency; a is the radio spectral index, 
which is connected with the value of the spectral index of 
the particle distribution, q, as a = {q — 3)/2, a(a) is the 
coefHcient fully depending on the value of a (Ginzburg 
1981; p.93). Assuming the magnetic field is dependent on 
the density as i? ex p*^, we represent its radial distribution 
inside the remnant as 



B = Bo-[x{r,Ms)Y 



(19) 



where Bq is the ISM magnetic field strength, and x{r, Ms) 
the radial density distribution (A. 3). 

Taking the definitions (6)-(8) into account, the 
badly determined parameter of poi can be expressed as 
(P0ic)2 = (V'PToC)2 = ^2 . 2mec2To«, with Tos = 3.13 X 
10^^/iVQ(j (org), where voo = vos/(10'' cm/s) is the initial 
shock velocity when the shock radius was Rq. After sub- 
stituting (17) into (18) and taking the last considerations 
concerning B and po into account, we obtain the expres- 
sion for the emission coefficient in the form 



where 

= 9.97 X 10-3^ r?_3 n^o B_5 2^*= 4^ 
(erg cm~^sr~^Hz~^), 



(20) 



(21) 



and 



_ a{a,) q, [:r(r)]^("'+i) 
^ ~ a(0.5) 4 ^ 23fe 



X (3.21- 10"^ ci)"' 



ai-0.5 / Mi 



(22) 



where Ci = fitp"^ Vqq B-^/pq , rj-3 is t] in units of 10~^, i?_5 
is -Bo in units of 10~^ G and ug is the radio frequency 
in units of lO'' Hz. The radially dependent part of the 
emission coefficient (22) is constructed in such a way that 
at r = 1 and Ms = Mo it becomes unity. It is important 
to note that the poorly known parameter ijj enters the 
formula for emissivity through ci as a combination with 
other parameters and, as shown below, the dependence of 
the evolutionary part of the emission coefficient on Ci is 
very weak. 

Using (20)- (22), the formula for the radio flux density 
can be written as 



Sv — Sqv ■ /i/(i?s) 

where 

So. = 387 r/_3 rieo 2^'= B_5 • -Ropc c'kpc (Jy) 



(23) 



(24) 



3 ^ 

and the surface brightness 
where 



So. = 1.23 X 10-i«r;_3neo 2^'= S_5 cj-^ii, 



(W 



Opc 

m-^sr-^Hz-i) 



and 



r dr. 



(25) 



(26) 



(27) 



(28) 



In (24) and (27), i?opc is the initial radius of SNR in pc 
and dkpc is the distance to SNR in kpc. The ratio Rs/Ro in 
(25) and (28) can be expressed through the Mach number 
by using the approximated formula (A. 2). 

2.3. Common properties of the model 

In this section we study general properties of the model, 
namely, the evolution of the main radio characteristics 
(emissivity, luminosity, surface brightness, spectral index) 
of the SNR with Mach number and their dependence on 
the input parameters of the model. The Mach number is 
used as the parameter describing the evolution of the SNR 
although it can bo easily transformed into other measures 
of evolution such as radius and/or time by using Eqs (A. 2) 
and (A. 6). There are three input model parameters: initial 
Mach number Mo, the power k in the relationship B oc p'' , 
and the combination of parameters Ci = p, V'^Vgg • B^^/vg. 
Both, Mo and Ci are not fully independent, but in this 
subsection we adopt them as independent parameters. The 
value of ci is restricted within the range between « 1 and 
pa 200 when the variables vgg and -B_5 are changing in the 
range 0.5 — 5 and 0.3 — 1, respectively. It is important to 
note that, as shown below, the main radio properties of 
the SNR depend very weakly on the value of ci . 

In Fig. la the normalized radial profiles of the volume 
emissivity coefficient, ei,{r), are shown for different sets of 
input parameters. As seen in this figure, emission at high 
Mach numbers is highly concentrated directly behind the 
shock front; but with decreasing the shock strength, the 
profiles broaden and, ultimately, at some critical Mach 
number, Ma, the peak of the profile separates from the 
shock front. This occurs around A^cr ^ 5 when the de- 
crease in the efficiency of acceleration due to steeper spec- 
trum of the newly accelerated electrons becomes more ef- 
fective than the adiabatic cooling of early accelerated elec- 
trons. As the magnetic field is frozen well in the ionized 
plasma, its strength reaches its maximum directly behind 
the shock front, but the peak in the distribution of radio 
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Fig. 1. The evolution of the radial profiles of an undimensioni 

coefficient of ciiiission £i/ (r) (Eq.22) in its dependence on tl: 
input parameters. 

emitting electrons is displaced from the regions where the 
magnetic field strength is at its maximum. Thus, the ef- 
fectiveness of transforming of the energy of particles and 
magnetic fields into the synchrotron emission decreases 
(Asvarov & Guscinov 1991). 

The effect of broadening the profile is more prominent 
for high values of A^o and small ci at fixed other param- 
eters. The radial distributions of emissivity presented in 
Fig. 1 lead to the radial distributions of the radio inten- 
sity on the face of SNR shown in Fig. 2 for several set of 
input parameters calculated by integration along the line 
of sight I{z) = J £j/(r) • dl, where z is the radial position 
on the radio image of SNR. 

The main features that can bo seen in this figure are: 

1) In evolved SNRs the radius of the radio remnant de- 
fined as the radius of the circle of the maximum inten- 
sity is smaller than its real radius, though for young rem- 
nants the peak of the intensity occurs directly behind the 
shock front. This implies that for old adiabatic SNRs with 
A^s < 10 the radio diameter systematically smaller then 
its real diameter. The difi^erence between these diameters 
increases with decreasing shock Mach number. 

2) As can be seen in Fig. 2, broadening of the radio shell 
in the course of the evolution of the SNR is the next prop- 
erty of the model. 
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Fig. 2. Normalized projected radial profiles of synchrotron 
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Fig. 3. The dependence of undimensional luminosity (Eq. 
25) on Ma. The tracks are labeled with input parameters 
{Ma, k, ci) 

3) The center-to-limb ratio of the projected profile of ra- 
dio emission is the measurable characteristic of an SNR. 
Our model predicts for this ratio a value of ^ 0.4 at high 
Mach numbers and increases up to a value of ~ 0.5 with 
decreasing Ms- At given Ms and Mq the dependence of 
the shape of I{z)— profiles on the other input parameters 
is weak. 

Figure 3 shows the evolution of the undimensional lu- 
minosity C with Ms for several sets of input parameters 
(A^o; k; ci). This quantity at first increases as li, oc M~^ 
with decreasing Mach number; then after reaching its 
maximum at Mmax, it decreases as oc Mg'^ when 
fc = 2/3 and as ~ At g for A; = 1. The corresponding 
dependencies of l^, on the shock radius Rg also can be 
approximated by two power-law relationships: l^, oc R^-^ 
and t oc R-^-^{k = 2/3) or ~ R-^-'^{k = 1). The values of 
Almax and the corresponding peak values of the undimen- 
sional luminosity Ivmax. are given in Table 1 . The value of 
^i/max linearly increases with Mo but the dependence on 
the parameter ci is weak. 
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Table 1. The parameters of the model, amax is (Jv at Ain 
uio is av at Ma = 10 
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1000 
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7.9 
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500 
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100 
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100 
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Fig. 4. The dependence of undimensional surface brightness 
(Eq. 28) on Mb- Curves are labeled with input parameters 
(A^o, k, Cl) 

Figure 4 shows the evolution of undimensional sur- 
face brightness, a^, (see Eq. 28) with A^s for several 
set of input parameters. As follows from the analysis of 
tracks, the relationship Gi, — M.s has a universal character 
within the framework of the adopted model. The depen- 
dence a I, = (7,y(A4s) can also be approximated by the two 
power-law dependencies: o"^ oc A4g '^^°-* at M^ > 10 and 
a^, oc A^s~^ at M.s < 5 — 6. The transition to declining 
branch is smooth and begins at Ms < 10. The values of 

at Ms = 10 and Ms = Almax are given in Table 1. At 
low values of A^s, the slope of the curve = (Jv{M.s) is 
steeper when the degree of relationship between the mag- 
netic field and the density is stronger (curve with fc = 1 in 
Fig. 4). At high Ms^ the value of depends weakly on 
Mq and k and barely depends on parameter c\. The ratio 
of the extreme values of calculated for the set of input 
parameters within their maximum plausible boundaries is 
a factor of 2.5 at Ms = 10 and does not exceed a factor 
of 3.5 for A/Js > 4, whereas this ratio for 1^ is much higher 
(see. Tabic 1). 

Although this is a very important prediction of the 
model, the Mach number is not an observable parameter. 
Remnants with the same values oi Ms may have com- 
pletely different diameters. Converting the dependence of 
the surface brightness on the Mach number into the sur- 
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Fig. 5. The evolution of the integral radio spectral index with 
shock Mach number. Curves are marked by the values of par 
rameters MQ,k,c\ 



face brightness - diameter dependencies, we find that at 
A^s ^ 10 surface brightness changes with SNR radius as 
j^-(o.3-o.5)^ followed by a steeper decrease with radius as 
Rs ''^'^ ' at Ms ^ 8. This transition between two power 
laws is smooth enough, because it begins at Ms ^ 10 
and finishes at Ms 5. Moreover, if we take into ac- 
count the above-mentioned effect of lagging of the radio 
shell behind the main shock front, the slope will be even 
steeper. The prediction of the model about the surface 
brightness-diameter relationship agrees with the conclu- 
sion of Berkhuijsen (1986) that SNRs evolve at a constant 
surface brightness followed by an abrupt drop, although 
our conclusion is true if the evolution of SNR takes place 
without the onset of the radiative phase. The case with 
possible radiative cooling is considered below. Our result 
disagrees with the conclusion of Berezhko & Volk (2004) 
that in the Sedov phase Sr cx Rs which is evidently 
the result of the assumption made by these authors about 
generation of the magnetic field at the shock front. 

The radio spectral index is the most important pa- 
rameter of SNRs. Our model gives very specific predic- 
tions for this parameter. The evolution of the integral 
spectral index with Ms is plotted in Fig. 5. The value 

of the spectral index is calculated using the formula: 

1 1 
a = J a{r) ■ £^{r)r'^dr / J ey{r)r'^dr for any given Ms- 

Tin Tin 

First of all, as can be seen in Fig. 5, the shape of all curves 
is very similar; the spectral index increases with decreas- 
ing Ms (reaching the value of ^ 0.6 at Ms < 4) until it 
achieves its maximum around A^g = 3 and then slowly de- 
creases. Although the shape of all curves are similar, very 
little depending on the input parameters, the maximum of 
the spectral index ctmax depends on the input parameters 
of the model, namely, ctmax is larger for small Mo at fixed 
other input parameters and it increases with increasing 
Cl. If we take the quadratic relationship between ci and 
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^Ao into account, the last conclusion implies that amax is 
higher for high values of the magnetic field strength. 

The general properties of the model that have prac- 
tical importance are the shapes of a^, = a^{Ais) and 
a = a{Ms) dependencies and the fact that they have a 
universal nature in the sense that they barely depend on 
the poorly known injection parameters. 

3. Application of the model to the real SNRs. 
T, — D relationship 

In the study of SNRs the empirical T, — D dependencies 
play an important role. This relation is used as the main 
tool for testing the results of the theory. Moreover, many 
statistical studies of Galactic SNRs have relied on the 
T, — D relations to derive distances for individual SNRs. 
The results of the previous section are presented in a func- 
tional form with arbitrary normalization and, also, with- 
out taking into account the possible onset of the radiative 
phase in the course of SNR evolution. In this section we 
apply these results to the real SNRs evolving in real in- 
terstellar environments. 

3.1. Input parameters 

To apply our model to the real SNRs, we first have to 
specify additional new parameters concerning the SNR it- 
self (the ejecta mass Mej, the initial kinetic energy of the 
ejecta Esn, the rms velocity of the ejecta vq), and the 
ISM (the total pressure Pq, the density of plasma rieo, the 
magnetic field strength Bq). 

The beginning of the model is parameterized by the 
initial radius Rq and Mach number AIq of the shock wave 
at which the Scdov phase begins. We ignore the ejecta 
stage because the possible contribution of this part of the 
remnant to the total emission with time will rapidly de- 
crease. 

In the following we assume a fully ionized gas with a 
ratio of specific heats 7 = 5/3 and with a helium abun- 
dance relative to hydrogen as riHe = 0.1 nn- The number 
density and electron density are n = nue + nu = 1.1 riH 
and rio = tt-h + Znjjo = 1-2 nn, respectively. The mean 
mass per particle and per electron in units of the proton 
mass, TOp, are fi = (1.4/2.3) and /io = (1.4/1.2), the total 
number of particles per electron is 2.3/1.2 « 1.92. The 
mass density of matter we express in terms of the electron 
concentration rie- p = PcmpUc. 

To take very extended SNRs into consideration, we use 
an analytical approximation of Cox & Anderson (1982). 
The details of this approximation are given in Appendix 
1. The pressure of ISM during the late stages becomes 
a leading factor in defining the evolution of the SNR, so 
it is important to include the pressure of the magnetic 
fields Pm = P^/Stt in the total pressure of matter in the 
ambient ISM and take the effect of the magnetic field into 
account in determining the value of the compression ratio 
at the shock front. For simplicity we make no difference 
between the thermal gas and cosmic ray components when 



we determine the total pressure in the ambient medium, 
and we assume that there is an equipartition between the 
three components of the ISM. 

Since the maximum speed of propagation of small 
perturbations in the magnetized medium is the speed of 
magneto-sound waves, therefore, as a measure of intensity 
of the shock wave, we use the ratio of the velocity of the 
shock wave Vg to the speed of fast magneto-sound waves 
Cms = (Cg + v^)""^/^ (where va is the Alfven speed, Cg the 
sound speed) in ISM : Ai = Vg/cus and for simplicity, we 
refer to this value as the shock Much number. Using the 
shock jump conditions for the perpendicular plane-parallel 
shock, we can find the compression ratio as the real root 
of the quadratic equation (Field et al. 1968): 



2(2-7) ^: 



2 + 7/3 



T7^+ 7-1 



27(/3+l) 1 
2 + 7/3 M"- 



(7 + 1) = 0, 



where f3 = Pth/Pu (Pih is the sum of the thermal gas and 
cosmic ray pressures) - one of the main input parameters 
characterizing the ISM in our model. This equation pro- 
vides a necessary condition that Xg ^ 1 with Ai ^ 1; 
and when /3 — + oo {M is ordinary sonic Mach number), 
it gives for the compression ratio the well-known relation 
(Landau & Lifshitz 1986): 



Xs{M) 



7+1 



7 - 1 + 2A4- 



The total pressure of the ISM can be expressed in terms of 
the interstellar magnetic field strength and the parameter 



Po4 = 2.88 (/? + 1)B% 



(29) 



where P04 is Pq in units of 10^ Kcm ^. For the initial 
Mach number we now obtain 



Mo = 458 



V5/3/6 + 1 



1/2 



V09 



(30) 



where vog is initial shock velocity in lO^cm/s. The mass 
of the material ejected at the SN explosion, Mej, and its 
kinetic energy, Ssn, are used as input parameters of the 
model characterizing the supernova event. The relation 
between them and vog is 



i?05l/v 



(31) 



where E051 is £^sn in 10 erg, and Wej is Mcj in solar 
masses. 

A very important parameter of the model is the ini- 
tial shock radius Pq at which the Sedov phase begins. 
Usually the beginning of the Sedov phase is determined 
from the condition that swept-up mass becomes equal to 
the ejected mass, e.g., at x = 1 in 4ttRqPq/3 = X'^ej- As 
the transition from the free-expansion phase to the Sedov 
phase is smooth and occurs over a finite period of time, it 
is difficult to choose the concrete value for Rq. Moreover, 
as shown by Cioffi et al. (1988) at some conditions (not 
far from the typical ones for real SNRs), the interior of 
the adiabatic blast wave may not settle completely to the 
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Sedov similarity solution. We determine the value of Rq 
from the law of motion (Al) and its new representation 

(A2): namely, equating them wc have (36/5)i?j? = RqA4q 
from which taking into account the definition of the ini- 

1/2 

tial Mach number as Ado = vq/ {iPq/Po) and (31), we 
obtain 



Ropc = 1.87 • (mej/Merieo)^/^ 



(32) 



whore i?opc is Ro in pc. The comparison of this formula 
with the one derived by introducing the parameter %, 
Ropc = 2.13 • {xiT^c}/ f^cn-coY^^ , gives x ~ 0.6. In our case 
the low values of x f^i't- physically acceptable because the 
self-similar structure in the shock interior begins to form 
earlier in the free expansion phase. Initially it occupies a 
small part of the volume behind the shock front, but with 
time the fraction of the volume with self-similar structure 
will increase. As mentioned above, we consider only the 
volume between Rin and Rs, where the Sedov structure 
holds, and ignore the emission from the inner volume of 
the remnant, which is occupied with ejecta material. 

When the downstream temperature approaches some 
"cooling temperature" Tc ~ 10^ K, the remnant will 
make a transition from the Sedov regime to the radiative 
or "pressure-driven snowplow" phase (McKee & Ostriker 
1977; Cioffi et al. 1988). A cold dense shell begins to form 
when the first portion of the shocked gas has radiated most 
of its thermal energy. Since we consider the acceleration 
of electrons from the high-energy tail of the downstream 
Maxwellian distribution function, it is natural to assume 
that injection of fresh particles into the process of accel- 
eration stops when radiative cooling of the shocked gas 
first starts. The velocity at the beginning of the radiative 
phase according to CiofR et al., (1988) is 



VPDS 



413 Eq^I^ nl^'^ km/s. 



The process of the radiative cooling was studied by many 

authors, and their results did not always agree. For in- 
stance, the same functional formula is derived by Franco 
et al. (1994), but with a lower value of the coefficient of 
265. In the following we employ the formula of Cioffi et 
al. (1988) because the authors of this work have given a 
more precise description of the preceding evolution of the 
shock wave and considered the very beginning of shell for- 
mation. In our notations the above expression takes the 
form 



(33) 



The shock wave Mach number at which thin shell for- 
mation begins is 

A^eooi = 18.42 (5/3/6 + 1)-^/' Me^' Elil%%'' Bzl (34) 

Prom this formula we can find the density of the ISM at 
which SNR will evolve without shell formation up to some 
critical value of Mcr as 



rieo ^ 1.08 X 10" 



7/9 pl/9 
A'e -^051 



(35) 



As the pressure in the ISM is the quantity that is to a 
lesser degree subjected to strong fluctuating, it is useful 
to rewrite this expression in terms of this quantity. Using 
Eq. 29 we have 



neo< 4.74X10 ^^^j 



(36) 



The SNR with E051 « 1 evolving in the ISM with 
/3 PS 2 (equipartition between thermal, CR and mag- 
netic pressures), /io ~ 1.4/1.2, and Po4 ~ 3 (Boulares 
& Cox 1990) will finish its life by merging with the 
ISM {Mcr ~ 2) before cooling becomes important, if 
rico < 2.65 X IQ-^ cm-^. As shown above, the break 
in the S — _D track to a steeper power law starts at 
A^break ^ 10 and, for this to happen without dense shell 
formation, the density of the ISM must satisfy the condi- 
tion rieo < 0.23 cm"^. The constraint on the density is not 
very strong, so, for many remnants in the course of their 
evolution the break in the T, — D track can occur prior to 
the beginning of the radiative phase. 

As a result, if the density of the ISM satisfies the value 
given in Eqs. (35) and (36), then the break in the S — 
D evolutionary track occurs due to the steeper spectrum 
of the newly accelerated electrons, which happens at a 
diameter 



i'break ~ 80.9 (EqSi/ ^04)'^^ (pc). 



(37) 



In the opposite case, when SNR evolves in higher density 
ISM, the brake in the T, — D track occurs due to the onset 
of the radiative cooling at 



i'break « 36.4 (x/Me)'/' E^l uj'^ (pc). 



(38) 



In deriving Erj^s. 37 and 38, we have assumed Albroak = 10, 
(3 = 2. Note that in the case of radiative cooling, the track 
has sharp break to steeper power law but, in the case of 
DSA, the track steepens gradually, and it is difficult to 
determine the precise diameter fbreak and Eq. 37 deter- 
mines the diameter of the remnant at which the surface 
brightness just begins to break. 

Finally, two additional parameters concerning the 
mechanism of acceleration r/ and remain as free parame- 
ters of the model. In principle, special radio, X-, and 7-ray 
observations of an individual remnant can be used to es- 
timate these parameters. Here we aim to examine the ca- 
pability of DSA to explain the averaged radio properties 
of shell-type SNRs and to reveal the principal manifes- 
tations of the action of DSA in SNRs. As shown in the 
following section, the statistical S — _D relation can be 
used to obtain an order-of-magnitude estimation for these 
parameters. In model calculations for real SNRs, we have 
used an estimation of 77 = 6 x 10"^, which is derived there. 

3.2. The T, — D relationship 

Empirical T, — D relations are very useful tools for test- 
ing various theoretical models of SNR evolution. A brief 



Asvarov A.I.: Radio emission from SNRs 



11 



review of this problem is given in the recent work by 
Urosevich et al. (2005). In Fig. 6 the surface brightness at 
the 1 GHz(EiGHz) ^ diameter diagram is presented, which 
consists of 158 large-diameter SNRs and candidates for 
SNRs with known distances: 33 composite and shell- type 
Galactic SNRs, 4 Galactic Loops, 25 in LMC, 7 in SMC, 
51 in M33, 30 in M31, 3 in NGC 300, 2 in NGC 6946, 2 
in NGC 7793, and one SNR candidate in IC 1613. For the 
galactic SNRs, we used only the data given in the cata- 
logue of Green (2004a, 2004b) with one exception - we did 
not include the SNR OA184 in our consideration because, 
as shown in the recent work by Foster et al. (2006), this 
object is most probably an HII region, not an SNR. The 
data on the Galactic loops are from Berkhuijsen (1986). 
The data for the remaining 121 SNRs in nearby normal 
galaxies are from the work of Urosevich et al. (2005). We 
did not use the remnants in the stur burst galaxies because, 
on the one hand, the majority of them arc very young ones 
and, on the other, the physical conditions in the ISM of 
these galaxies differ from the standard conditions in the 
normal galaxies considered in our model. 

In the sample of SNRs presented in Fig. 6 we also did 
not include Radio supernovae and plerionic SNRs because 
the mechanism of their emission differs radically from that 
of our model. Our model also cannot explain the radio 
emission from the young SNRs actively interacting with 
the dense circumstellar structures or else in the double- 
shock stage of development. For instance, the radio emis- 
sion from the Cas A SNR is commonly considered as the 
result of the interaction of the SN ejecta with circumstel- 
lar material. The radio emission from other young histor- 
ical SNRs Tycho, Kepler, and AD1006, in principle, can 
be explained in the framework of our model. Further, our 
model meets with difficulty in explaining the radio emis- 
sion from SNRs with power-law spectra that are too flat 
or too steep. 

Some evolutionary tracks are also shown in Fig. 6 and 
the parameters of the models used to calculate them are 
given in Table 2. In all runs, we used an estimate of 
6 X 10^* for the injection parameter rj, derived as follows. 
For each value of surface brightness in the Siqhz — D dia- 
gram, the largest observed diameter is considered to be the 
value of fbreak duc to the radiative coohng (Eq. 38). The 
largest diameter remnants for a given surface brightness 
occupy the upper right region in the empirical Siqhz — D 
distribution, which should be reasonably free of selec- 
tion effects. For reliable estimates, it is desirable to con- 
sider remnants with intermediate values of E because the 
brightest, largest SNRs can be the ones for which our 
mechanism is not applicable and, in the opposite case, 
the largest SNRs with low E are affected by selection ef- 
fects. As seen in Fig. 6, Galactic SNR Kes 67 (G18.8-h0.3) 
with IgEiGHz = —19.58 and lgl?(pc) — 1.75 satisfies the 
selection requirements. As a turning point of the track, 
along which this remnant evolves, we take the point with 
coordinates IgEicHz = —19.55 and IgD(pc) = 1.75 be- 
cause, besides Kes 67, one more remnant, B0547-697 in 
LMC with IgEiGHz = -19.51 and IgD(pc) = 1.75 , is 



Table 2. Input parameters of the model runs. iJosi is -Esn in 
units of 10®^ erg, B_6 is So in units of 10^* G, Dbreak is in pc. 



all rans 


k = 


2/3, »7 = 


6 X 10- 


-^/3=2, 


mej = 


1.4 


No 






B-5 


Dbreak 


Mcool 


Cl 


1 


5.0 


1 


0.6 


22 


37 


5.40 


2 


1.0 


1 


0.6 


46 


13 


5.40 


3 


.5 


1 


0.5 


60 


10 


4.50 


4 


.1 


1 


0.5 


63 


3 


4.50 


5 


.05 


1 


0.3 


88 


4 


2.70 


6 


.005 


1 


0.3 


89 


1 


2.70 


7 


.001 


1 


0.3 


88 


1 


2.70 


8 


.5 


0.1 


0.5 


29 


9 


0.45 


9 


.1 


0.1 


0.5 


29 


3 


0.45 


10 


.01 


50 


0.3 


323 


2 


135.0 



located around this point. Using these data in Eq. 38 we 
find rieo = .32 cm~'^. Then, using the latter value in Eq. 27 
and adopting B-5 w 0.5 and £^051 = 1, we get an estima- 
tion of ?7 « 7.0 X 10~^. Varying the values _B_5 and £Jo5i 
within reasonable limits 0.3-^ 1.0 and 0.5 -j- 2, respectively, 
we find a range for 77 of (3-^11) x 10~*. In the calculations 
we have used 77 = 6 x 10""*. 

First of all, it is important to note that the shapes of 
the model's evolutionary tracks are in very good accor- 
dance with the claim made by Berkhuijsen (1986) that 
SNRs may evolve at nearly constant E followed by a steep 
decrease. In our model such behavior of the E — £> tracks 
has a natural explanation. In the low-density ISM, when 
SNR evolves without the onset of the radiative cooling, 
this is the direct consequence of the action of DSA, but 
when SNR c;volvc>s in the high density ISM, it follows from 
the assumption that acceleration of electrons takes place 
from the downstream thermal distribution. This assump- 
tion also naturally results in the tight correlation between 
the surface brightnesses in radio and X-ray wavelengths 
revealed by Berkhuijsen (1986). 

From the shapes of the evolutionary tracks and the 
very wide possible range of variation in the surface bright- 
ness for a given diameter (up to 3 orders of magnitude), 
it follows that there is no physically meaningful Y, — D 
relation that can be used as the tool for estimating the 
distances to SNRs. For example, two identical SNRs evolv- 
ing in identical environments but with different diameters 
(and ages) may have very close values of E. On the other 
hand, for two SNRs with the same D the values of E may 
differ by several orders of magnitude. 

To see this, let us write the expression for E (Eq. 27) 

as 

where it is assumed that in the ISM the magnetic field is 
correlated with the density as Bq oc tIq" , which can easily 
give a range of variation in E up to 3 orders of magnitude. 
Note that in the work of Berezhko & Volk (2004) , the value 
of Er at a given diameter depends only on the SN energy 
as £/gN ■ 
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Fig. 6. EiGHz — diagram for 158 com- 
posite and shell-type SNRs and SNR 
candidates in our Galaxy (filled circles) , 
M31 (small crosses), M33 (open dia- 
monds), LMC (open triangles), SMC 
(filled triangles), NGC 300 (filled dia- 
mond), NGC 6946 (x 's), NGC 7793 
(open circles), and IC 1613 (aster- 
isk). Galactic loops are shown as large 
crosses and the only candidate for HNR 
as a large open circle. For several galac- 
tic SNRs, error bars due to the distance 
uncertainty are shown, and for one 
galactic SNR (HB 9) only upper limit 
for the distance is known. Theoretical 
S — D tracks are labeled according to 
the numeration in Table 2. 



Concerning the "maximum observable diameter" dis- 
cussed by Berkhuijsen (1986), our model gives the relation 

EiGHzOci?^^/^«+^°D~iW'=o/2); (40) 

in particular, at fco ~ 0.5 we have the relation E cx -D,7iax^ 
derived in (Berkhuijsen 1986). 

The observational 'E — D relations, usually presented in 
the power-law form (E oc D^), are mainly due to: 1) pres- 
ence of SNRs with diameters from a wide range of sizes 
from small-diameter bright young SNRs and large faint old 
SNRs; 2) effect; 3) selection effects. Unfortunately, it 
is difficult to check the influence of various selection effects 
on this problem. For instance, the absence of a statisti- 
cally meaningful E — _D relation for the sample of SNRs 
in M33, which is particularly incomplete in the smallest 
and extended faint remnants, was emphasized by Gordon 
et al. (1999). In Fig. 7 the 1.4 GHz flux density - diam- 
eter relation is presented for 121 SNRs in nearby normal 
galaxies if they are placed at 1 kpc distance. All the data 
are from Urosevich et al. (2005). For an equidistant sam- 
ple of remnants, this relation can be used as an equivalent 
to the E — D relation, but it is free of the D^^ effect. As 
noted above (see Sect. 2.3) the effect of D^^ results in a 
smaller difference between the different Yi — D tracks than 
it is for luminosity - diameter tracks with the same set of 
the input parameters. 

Thus, in the formation of E — relation, the two 
main parameters no and i?sN have equal importance (See 
Eq. 39). At intermediate diameters, the contribution of 
no is more significant, but £^sn critically determines the 
maximum size of the radio remnant, especially, in low- 
density environments. As can be seen from Fig. 6, practi- 
cally all the remnants can be described within the range 
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Fig. 7. 1.4 GHz flux density - diameter relation for 121 SNRs in 
nearby normal galaxies -I- Cas A SNR, if they are at 1 kpc dis- 
tance. Symbols and lines have the same meanings as in Fig. 6. 

of (10^*^ 10^^) erg. Varying mainly the value of E'en, 
the model is able to explain the existence of very large 
diameter SNRs and even the remnants of extremely ener- 
getic SNe, known as hypernovae (HNRs). For instance, our 
model easily explains the radio emission from the Galactic 
loops, so there is no need for energies higher than the stan- 
dard value of ~ 10^^ erg. On the contrary, according to 
the model for the existence of such a large radio source as 
the SNR candidate N7793-S26 in the galaxy NGC7793, at 
least '-^ 5 X 10^^ erg energy is required. Although in Fig. 6 
only one track passing through this object is shown, of 
course, this is not the only possible curve; but for any 
plausible set of other input parameters, this constraint on 
energy is indispensable for explaining the radio emission 
from this object. Such energies are characteristic of hy- 
pothetical Hypernova explosions, and this result can be 
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considered as evidence of the real existence of such ob- 
jects, and the SNR candidate N7793-S26 can be regarded 
as a real HNR. It is important to note that no evidence 
has been found of an interior star cluster, the detection of 
which would argue in favor of the collective SNR hypoth- 
esis (Pannuti et al. 2002). The idea that the detection of 
the very large diameter radio remnants may be considered 
as proof of the real existence of hypernova was proposed 
in our early work (Asvarov 2001), where the parameters 
of expected hypernovae remnants were calculated. As seen 
in Fig. 6 at intermediate diameters lOOpc), the adia- 
batic radio remnants of Hypernovae, evolving in the ISM 
with typical parameters, do not differ in essential features 
from the ordinary SNRs. Of course, in the early stages of 
their evolution the hypernovae can be very bright in radio, 
depending on the properties of the circumstellar matter, 
but at the same time, the high luminosity of a very young 
SNR does not mean that we are necessarily dealing with 
a hypernova. 

We must note that in terms of other known models, 
e.g., the model of van der Laan (1962) and its modifica- 
tions, describing the radio evolution of evolved shell-type 
SNRs, it is not easy to explain the existence of SNRs with 
such large diameters 

Within the framework of the present model, a number 
of features of the empirical T, — D relation has a simple 
explanation. For instance, the small number of remnants 
with small diameters and low S (lower left corner in the 
diagram) is the result of the very fast evolution of the 
SN blast wave in the low-density ISM where SNRs have 
low E. According to our model, it is easy to account for 
a high concentration of SNRs at diameters of 20 — 50 pc 
in the S — Z? diagram: various evolutionary tracks inter- 
sect at these diameters and the sample of remnants here 
consists of objects evolving at different initial values of 
-BsN and noe- However, it seems more likely that this ex- 
cess comes from the energies of most SNRs being closer to 
Iq50 ^j.g ^}jan 10^^ erg (Lozinskaya 1980). Another fact in 
favor of such a possibility is the absence of SNRs near the 
break points of 10"''^ erg tracks in a high density medium 
(curves (1) - (3) in Fig. 6), although this part of the E — D 
distribution is free of the selection effects. Unfortunately, 
the importance of selection effects at large diameters is 
hardly known, but if most SNRs have initial energies of 
~ 10^^ erg, then it follows from the comparison of the 
model with the E — I? distribution that there should be 
a large population of 150 — 250 pc diameter SNRs with 
SiGHz < 10"^^ Wm~^ sr~^ Hz~^. Identification of such 
low-brightness remnants in our Galaxy is a very difficult 
task, although there are a few successful attempts (Reich 
2002); but they can be detected in nearby galaxies. 

From the distribution of the model tracks in the S — D 
diagram, a very important conclusion can be made that 
shell-type SNRs with Sighz < IQ-^o Wm-^sr-^Hz-^ 
are located in the low-density phase of the ISM. It is 
important to take into account the fact that our sample 
of Galactic SNRs consists of remnants for which the dis- 
tances are mainly determined from the interaction of the 



supernova blast wave with the dense environments; i.e. 
our sample consists of objects located in special environ- 
ments, e.g., near large atomic or molecular clouds. Indeed, 
for 40 Galactic SNRs with estimated distances out of 101 
SNRs (without plerions), for which the Eiqhz can be cal- 
culated in the catalogue of Green (2004a) , the mean value 
of the surface brightnesses Eiqhz (in standard units) is 
4.82 X 10-19 ( 7.47 x IQ-^o without Gas A), but the mean 
value of El GHz for remaining 61 SNRs (without distance 
estimations) is 7.45 x 10^^^. 

Finally, the fact that 63 (62%) out of those 101 SNRs 
have EiGHz < 10"^" Wm"^ sr^^ Hz^^ can be considered 
as a proof that the large fraction of SNRs are evolving in 
the low-density phase of the ISM. The same conclusion 
about the environments of the galactic SNRs was earlier 
made by Marsden et al. (2001) from the association of 
AXP and SGR with SNRs and in Koo and Kang (2004) 
from the considerable deficit of Hl-emitting, old SNRs. 
This result is also in accordance with the theory of McKee 
and Ostriker (1977), which predicts the widespread exis- 
tence of hot and low-density gas in the Galaxy. 

3.3. Spectral index 

The statistics of spectral indices is another argument fa- 
voring the DSA mechanism. The catalogue of SNRs of 
Green (2004a) contains 231 SNRs from which 180 (78%) 
are shell- type (labeled in the catalogue "s" and "s?") for 
91 of which the values of the spectral indices are deter- 
mined: for 87 SNRs the mean value of the spectral index 
is a = 0.506 ± 0.111 and 4 SNRs have variable spectral 
indices. Twenty-one of these 91 SNRs (23%) (7 SNRs with 
a > 0.65, 10 SNRs with a < 0.4 and 4 SNRs with variable 
a) have spectral indices that are not in agreement with 
DSA: our model predicts a in the range of (0.5 - 0.6). The 
occurrence of shell-type SNRs with a flat spectrum can 
be explained by the inclusion of second-order Fermi ac- 
celeration (Schlickeiser & Furst 1989; Ostrowski 1999) or 
by the effect of compression of the interstellar cosmic ray 
electrons by the radiative shocks (Cox et al. 1999). Indeed, 
there is some evidence that flat spectra are characteristic 
of SNRs interacting with higher density environments. 

Steepcir spectrum SNRs are as a rule the young ones 
(for example, the historical shell remnants Gas A, Tycho, 
Kapler, SN 1006 all have spectral indices steeper than 
a = 0.6), which can be explained by the nonlinear modi- 
flcation of the shock structure (see, for the recent review, 
Volk 2006), but in the case of evolved remnants there is 
no more or less satisfactory explanation for such a. In 
the framework of the DSA mechanism in the test parti- 
cle approximation, the values of a steeper than 0.6 can 
be explained if we assume that the remnant evolves in a 
pre-existing cavity with a radially increasing density dis- 
tribution. 

Our model does not predict any detectable dependence 
of the spectral index on the diameter because the univer- 
sal character of the a(A1s)-dependence (Fig. 5) does not 
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lead to the same character for the a{D) - dependence. 
But, in any case, the model predicts spectral indices ap- 
proaching the value of ^ 0.6 for large diameter SNRs. 
According to our model, the SNRs with Mach numbers 
Ms < 5 will have the mean value of a w 0.55 0.60. 
Assuming for simplicity that SNR evolves according to 
the Sedov law, Ms oc t~^/^, the number of SNRs with 
Mach number smaller than Ms is A''(< Ms) oc M7^^^ , 
from which it follows that the number of SNRs with 
a > 0.55, N{a > 0.55) = N{Ms < 5) must be about 
70% of the total number of SNRs. Here we have adopted 
M[ = 2.5 for the final Mach number at which S drops 
more than two orders of its initial value, and the SNR 
becomes invisible. Prom 92 SNRs only 18(20%) (exclud- 
ing 3 historical young SNRs: SN1006, Kepler, Tycho) have 
0.55 < a < 0.65. The shortage of SNRs with such spectral 
indices can be easily explained by various selection effects, 
all of which are biased against the faint large-diameter 
SNRs. 

3.4. Radiative SNRs 

Since electrons, by assumption, are injected from the high- 
energy tail of the downstream Maxwellian distribution 
function, the process of generating new radio-emitting 
electrons stops with the beginning of the radiative cool- 
ing. This implies that, with the beginning of the radiative 
phase, the radio emission from SNR starts to drop rapidly. 
On the other hand, with the onset of radiative cooling the 
emission from an SNR in X-rays will also decrease rapidly. 
Indeed, a tight correlation between radio and X-ray sur- 
face brightnesses has been revealed by Berkhuijsen (1986), 
which confirms our suggestion about the injection of elec- 
trons into the DSA. 

As a consequence of this, our model predicts the ab- 
sence of radio emission from the radiative remnants evolv- 
ing in the homogeneous ISM according to the standard 
theory of SNR evolution and the; existence of radio-quiet 
radiative SNRs emitting in optical and Hi wavelengths. 
The number of such objects is expected to be very large 
because the evolution of the remnants in the radiative 
stage takes place more slowly. Indeed, very recently first 
evidence of the existence of these SNRs was reported 
(Mavromatakis et al. 2005; Koo et al. 2006). As noted 
above, there is large deficit of radiative radio remnants in 
our Galaxy, and several optically bright SNRs are those 
interacting with the large-scale clouds, and our model does 
not describe such remnants. The deficit of radiative SNRs 
in our galaxy can be explained by the irregular high- 
intensity background emission, which complicates the de- 
tection of extended and often fragmented and clumpy 
shells of SNRs. On the contrary, in the case of extragalac- 
tic observations, much more optical SNRs have been de- 
tected (Matonick & Fesen 1997, Gordon et al. 1999). In 
this connection, it is important to note the very interesting 
finding by Duric and collaborators (Duric 2000; Pannuti 
et al. 2000; Pannuti et al. 2002) that very small amounts of 



optically-identified extragalactic SNRs are detected in ra- 
dio and X-ray wavelengths; and vice versa, among the ra- 
dio and X-ray selected SNRs, there are very few optically 
detected SNRs. The large number of extragalactic SNRs 
detected in optical wavelengths are partly the result of the 
fact that the sensitivity and resolution limitations severely 
reduce the effectiveness of radio and X-ray searches for 
many distant extragalactic sources. At the same time, we 
must also take into account that many SNRs located in 
and near HII regions easily can be missed in optical wave- 
lengths. Since the strong optical emission from an SNR or- 
dinarily implies that it is in a radiative evolutionary state 
(Ha is also observed from non-radiative remnants, such as 
Tycho's remnant, but the emission is weak), the absence 
of radio emission from the optically identified extragalac- 
tic SNRs can be considered a serious argument against 
the mechanism of van der Laan to be the main generator 
of nonthermal radio emission from the radiative SNRs. It 
thus seems very likely that most of the optically identified 
extragalactic SNRs represent an assembly of old remnants 
in a radiative phase evolving mainly in the dense ISM that 
have already faded or become weaker as radio sources, so 
these observational data confirm the conclusions of our 
model. 

An opposite explanation has been suggested by 
Pannuti et al. (2000) and Pannuti et al. (2002) that the 

optical remnants not possessing the radio counterpart are 
those evolving in the medium with low density. If this 
were the case, then these remnants would have consid- 
erably larger diameters compared with those that only 
radiate in the radio range. This is because the remnant 
should process larger volume to accumulate a column den- 
sity sufficient for the onset of the radiative cooling in the 
low density ISM . But the statistics of SNRs in M33 does 
not seem to support this possibility: the values of mean 
diameter of SNRs detected in radio and the optical SNRs 
without radio counterpart are almost equal: 37.34 pc (for 
44 radio SNRs from Gordon et al. 1999) and 38.81 pc 
(for 31 optical SNRs without radio from Gordon et al. 
1998), respectively. Unfortunately, it is difficult to take 
into consideration the influence of various selection effects 
on the statistics of extragalactic SNRs, but it is clear that 
many weak radio SNRs remain undetected. For example, 
the faintest radio SNR in the sample of SNRs in M33 has 
SiGHz = 8.4 X 10^^^ Wm^^ sr^^ Hz^^, which is an order 
of magnitude brighter than the faintest SNR in the sam- 
ple of Galactic SNRs, and from 101 Galactic shell-type 
and composite SNRs with known Eiqhz 11(11%) SNRs 
are fainter than that remnant in M33. 

4. Conclusions 

We have presented the model of the radio emission of shell- 
type SNRs evolving in the homogeneous ISM based on the 
assumption that the radio-emitting electrons are acceler- 
ated by DSA mechanism from the downstream thermal 
Maxwellian distribution. The aim of the work was not the 
detailed modeling of a certain SNR but to examine the 
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capability of DSA to reproduce the general statistics of 
shell-type SNRs and, as far as possible, to obtain con- 
straints on the mechanism for accelerating electrons. 

1. The main properties of the model that intrinsically 
reflect the action of the mechanism of DSA are as follows: 

— SNRs evolve at nearly constant radio surface bright- 
ness (S oc £)-(0-3+o.5)-j followed by a steep decrease 
(steeper than oc D~'^-^). The break to a steeper power 
law in the T, — D track begins at shock Mach numbers 
A^s ^ 10. This result is in excellent agreement with 
the conclusions of Berkhuijsen (1986). 

— With increasing diameter of the remnant, the radio 
shell broadens, and the center-to- limb ratio of the radio 
intensity over the surface of the SNR diminishes. 

— The radio size of the SNR with Adg < 10 is smaller 
than the real diameter of the shock front, and this 
difference increases with time. 

— The value of the integral spectral index of the remnant 
increases from the canonical value of 0.5 up to 0.6. The 
values of a for the evolved remnants that are higher 
than 0.6 can be explained in the framework of DSA 
in the adopted approximation, if the remnant evolves 
in a preexisting cavity with a density distribution in- 
creasing radially from the center of explosion. 

These features have a universal nature; i.e. they weakly 
depend on the input parameters of the model and prac- 
tically do not depend on the poorly-known injection pa- 
rameters. 

2. Our model explains many properties of the empiri- 
cal E — distribution for shell-type SNRs including the 
very large SNR candidates such as Galactic loops and 
N7793-S26 in the galaxy NGC 7793, which is considered 
as a real candidate for the hypernova radio remnant. From 
the shapes of the evolutionary tracks and the very wide 
possible range in variation of the surface brightness for a 
given diameter (up to 3 orders of magnitude according to 
Eq. (39)), it follows that there is no any statistically (also, 
physically) meaningful E — £> relation that could be used 
as the tool for estimate the distances to SNRs. 

3. Comparison of the model results with the observed 
E — _D distribution shows that most of catalogued galac- 
tic SNRs evolve in a low-density phase of the ISM which, 
in turn, may imply that the interstellar space in the in- 
ner Galaxy is largely filled with a very tenuous gas as in 
the three-phase ISM model. If this is the case, according 
to the model, a population of 150 — 250 pc SNRs with 
5^iGHz ^ 10~^^ Wm~^ sr~^ Hz~^ are expected to exist if 
the kinetic energy of the explosion is ~ 10^^ erg. However, 
the comparison of theoretical E— D tracks with the empiri- 
cal distribution shows that large fraction of observed SNRs 



have energies closer to ~ 10^" erg than to ~ 10^^ erg. 

4. The model considered here predicts that, with the 
beginning of cooling phase, the radio emission from the 
remnant starts to drop rapidly, i.e. we predict the absence 
of considerable radio emission from SNRs in the radia- 
tive stage of evolution. We interpreted this effect and sev- 
eral other observational facts as serious evidence that the 



mechanism of van der Laan is not an effective generator 
of the radio emission of radiative SNRs. Observational de- 
tection of a large number of radio-quiet SNRs would serve 
as confirmation of the model considered in this paper. 

5. By using the constructed empirical T, — D distri- 
bution, we obtain an estimation of the fraction of elec- 
trons accelerated from the thermal pool in the range 
(3-^11) X 10~*. If acceleration takes place directly from the 
high energy tail of the downstream Maxwellian distribu- 
tion function, then the corresponding injection momentum 
is estimated as pinj ^ (2.7 — 3) • pth • 

Appendix A: Cox-Anderson approximation 

The analytical approximation constructed by Cox & 
Anderson (1982) describes the development of an adi- 
abatic spherical blast wave in a homogeneous ambient 
medium of finite pressure. At early times, when the shock 
wave is strong, the structure is that of the usual Sedov 
self-similar solution. The law of expansion of the shock 
wave can be inferred from the expression (7 = 5/3) 



(3y. - 2)2 



2.02£'o 
75Po 



(A.1) 



which describes the evolution of the post-shock pressure 
P2 = Vs ■ Po with Rs] and at large t/s, it gives the Sedov 
solution: P2 = (3 • 2.02 £"0/25) R-^. Using the jump con- 
dition at shock front in the form yg = {5Mg — l)/4, this 
formula can be expressed in terms of Mach number as 



(M2 - 1)^ 



For our purposes it is convenient to use this expression in 
another form, namely, as 

R^.^R^.-Ml^^'-''/'f 



(A.2) 



which at Ais ^ 1 takes the form of Aig — 
M.o{Rs/ Ro)~^^^ , where Rq and Mq are the initial values 
of the shock radius and Mach number at the beginning of 

Sedov phase. 

In this approximation, the radial distribution of the 
density is 



.(r,E.).^ = 
Po 

X exp 



5 

2 ' 



■Q 



^9/2 
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(A.3) 



where = (4t/s + l)/(t/s -|- 4) = 4/(1 + ^/Ml) is the 
compression ratio at the shock front. The history of any 
element of matter can be evaluated by using the equation 



n(r,i?.)^|i=rV^ 



exp 



5/2 



(rO - 1) 



, (A.4) 



which connects the original location Ri of a mass element 
to its present position R = Rg ■ r when the shock radius 
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and Mach number are Rg and Ms, respectively. In (A. 3) 
and (A.4) 



2 x,(x,-l)(7x,-13) 

3 4x, - 1 

SM^M^s - l)(5.Ms - 13) 



(A.5) 



The expressions (A. 3)- (A.5) describe the structure of SNR 
completely. At A^s 1 they are in excellent agreement 
with the self-similar Scdov solution, and as noted by the 
authors, they were tested by the numerical calculations 
and are good up to j/g ~ 2 {Ms ~ 1.4). The dependence of 
Mach number and other dynamical variables on the time 
can be found through the equation 



1 = 1 

to 



{z"^- 1/5) dz 



3 ° (.2_l)2.(^2_n/15)V3' 



(A.6) 



which follows from Eq. (A. 2) and the definition Csdt = 
(1/A^s) (dRs/dMs) dMs- It is important to note that at 
high values of A^g this equation corresponds to depen- 
dence Ms = Mo[2.5{t/to) - 1.5]-3/^ but not Ms = 
Ma {t/to)~^^^ as in case of a pure Scdov law of expan- 
sion. Utilization of such a form for the dependence of Mach 
number on time corresponds to using the offset power-law 
formalism (CiofFi ct al. 1988), which describes the very 
beginning of the smooth transition from the phase of free 
expansion to Sedov phase. As a result, we can assign low 
values to the parameter x that allows us to cover a large 
part of the life of the remnant. 
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Appendix B: List of Symbols 

Bq - interstellar magnetic field. 

i?_5 - interstellar magnetic field in units of 10^^ G. 

C - parameter describing adaiabatic cooling of accelerated 

electrons [Eq. (12)]. 

ci = /xV'^Vgg • B-^/vg - parameter of the problem [below 
Eq. (22)]. 

Ds - SNR diameter. 

Dhreak " SNR diamctcr at vifhich break in the T, — D track 

occurs [Eqs. (37), (38)]. 

dkpc - the distacc to SNR in kpc. 

fa - distribution function of accelerated by DSA electrons. 

/ - radially dependent distribution function [Eq. (14)]. 

fm - MaxwcUian distribution function of electrons. 

H - Hcavisidc step function. 

-EsN - initial kinetic energy of the SN shell. 

£"051 - initial kinetic energy of the SN shell in 10^^ erg. 

fc, ko - power law index in the dependence of the magnetic 

field on the density B (x p'^ . 

Ke - coefficient of the power-law spectrum [Eq. (16)]. 
li, - undimensional radio flux density [Eqs. (23), (25)]. 
A4,A4s - shock wave Mach number. 

M cool ' Mach number at which thin shell formation 
begins [Eq. (35)]. 



Mr 



Mach number at which L reaches its maximum. 



Mo - initial shock wave Mach number. 
M, 



ejected mass at SN explosion. 
ruej - ejected mass at SN explosion in units of solar mass. 
Tie - electron number density, 
rieo - electron number density in the ISM. 
Po - total pressure of the ISM [Eq. (29)]. 
Po4 - total pressure Pq in units of 10'' Kcm~^. 
Pth - thermal gas -|- cosmic rays pressure. 



Pm - magnetic field pressure. 

P2 - post-shock pressure. 

Pinj - injection momentum of electrons. 

Pmax - maximum momentum of accelerated electrons. 

Pth - thermal momentum of electrons. 

PiO - injection momentum when shock radius was Rq. 

Q - parameter of the Cox-Anderson approximation [Eq. 

(A.5)]. 

q - spectral index of the momentum distribution function 

of accelerated by DSA electrons [Eq. (1)]. 

Qi - spectral index of the momentum distribution function 

of electrons accelerated when shock radius was Ri. 

R - radial coordinate from center of the SNR. 

Rc - characteristic radius of SNR [Eq. (A.l)]. 

Rs - SNR radius. 

Ro - initial radius of SNR at which Sedov phase begins. 
Ri - original radial position of mass element. 
Rin - present radius of mass element which was first 
shocked at i?o- 

— Rin I • 

ri - normalized original location of mass element = Ri/Rg 
[Eq. (A.4)] 

r - normalized radial coordinate from center of the SNR: 

r = R/Rs. 

Sv - radio flux density [Eq. (23)]. 
vos - initial shock velocity. 

vo9 - initial shock velocity in the units of 10^ cm/s. 

X - normalized radially dependent density distribution 

[Eq. (A.3)]. 

Xs - shock compression ratio ([Eq. 3] and Section 3.1). 
a - spectral index of the radio emission, determined as 

fi - ratio of thermal gas -I- cosmic ray pressures to the 
magnetic field pressure (Sect. 3.1). 

e,y, £1,0, £i'(r) - radio synchrotron emission coefficient 
[Eqs. (18), (20-22)]. 

e = Pinj/Pmax [Eq. (7)]. 

r) - injection parameter [Eq. (5)]. 
K - spatial diffusion coefficient. 

p, - mean mass per particle in units of the proton mass, 
/ic - mean mass per electron in units of the proton mass. 
S - radio surface brightness of the SNR. 
Soy - radio surface brightness of the SNR [Eq. (27)]. 
5^iGHz - radio surface brightness of the SNR at 1 GHz. 
(Tiy - undimensional surface brightness [Eqs. (26), (28)]. 
X - parameter characterizing the transition from free 
expansion to the Sedov phase. 
- injection parameter [Eq.(6)]. 



